MASSACHUSETTS INSTITUTE OF TECHNOLOGY 
ARTIFICIAL INTELLIGENCE LABORATORY 


A.I.Memo No. 1013 


February, 1988 


Utilizing Dynamic Stability 
To Orient Parts 


Neil C. Singer 
Warren P. Seering 


Abstract: The intent of this research is to study the dynamic behavior of a 

solid body resting on a moving surface. Results of the study are then used to pro¬ 
pose methods for controlling the orientation of parts in preparation for automatic 
assembly. Two dynamic models are discussed in detail. The first examines the 
impacts required to cause reorientation of a part. The second investigates the use 
of oscillatory motion to selectively reorient parts. This study demonstrates that 
the dynamic behaviors of solid bodies, under the conditions mentioned above, vary 
considerably with small changes in geometry or orientation. 
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1 Introduction and Problem Definition 


Conventionally, parts are oriented by bowl feeders. These machines vibrate and 
thereby convey parts through a series of filters which reject all but a particular orientation. 
Rejected parts are returned to storage. Many researchers have examined the implemen¬ 
tation of programmable or adjustable filter stages (Boothroyd, 1975) (Boothroyd and Ho, 
1977) (Boothroyd and Murch, 1970) (Boothroyd et al., 1982) (Boothroyd et al. 1977) 
(Lozano-Perez, 1986)(Murch, 1977) (Murch and Boothroyd, 1975) (Murch and Poli, 1977) 
(Redford et al., 1983a) (Redford et al., 1983b) (Singer, 1985) (Suzuki and Kohno, 1981). 
These techniques are often successful for limited classes of parts. The scope of this paper 
is to present theoretical results which may be useful in feeder designs. It is hoped that a 
more detailed understanding of the dynamic behavior of bodies resting on a moving surface 
will facilitate the design of new types of programmable parts feeders. 


2 Impact Reorientation 

In this section, we determine the conditions necessary to cause part reorientation. A 
part starts in a stable orientation (a natural resting aspect (Boothroyd and Ho, 1977)) on a 
fiat surface. It is then given an initial horizontal velocity after which it impacts a wall near 
the surface upon which it is resting (figure l). The equations describing the impact are 
based on conservation of angular momentum about the impact point. After impact, the 
principle of conservation of energy is applied in order to determine whether reorientation 
occurs. 

The part shown in figure 1 is a cylinder, but the analysis easily generalizes to any 
object. The important parameters are the mass, m of the part; the moment of inertia, 
/; and the constants, / and h, which locate the center of mass, cm, with respect to the 
impact point. The expressions for angular momentum, L,, before and, Lf, after impact 
are (Klepner and Kolenkow, 1973): 


Li = mV{h 

hVf ( 2 - 1 ) 

f ' 

where V,- is the speed of the center of mass prior to impact, Vf is its speed just after impact, 
and I a , the inertia about the impact point, is defined as 

I a — I c + {h 2 + l 2 )m 


with I c the inertia about the center of mass. Equating kinetic energy just after impact to 
the change in potential energy which will cause the part to reorient yields 
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Solving equations 2.1 for Vf, substituting into equation 2.2, and then solving the result 
for the critical value of V{ yields 


Vcrit 


\/2g ( h 2 + l 2 ) 2 - h [( h 2 + / 2 ) + 

h 


(2.3) 


which is an equation for the minimum initial horizontal velocity that will cause reorienta¬ 
tion. For cylinders, ^ is |h 2 + ~l 2 ; for rectangular parts, ^ is |(/i 2 + / 2 ). 

In order to test the ability of impact reorientation to discriminate among different 
orientations, it is useful to consider a part which has two very similar orientations. Figure 
2a shows a rectangular solid part in two possible orientations. One edge is bS longer 
than the other where is expressed as a fraction less than one. As the part exchanges 
orientations, the center of mass location changes in both height and horizontal distance 
from the impact point by 6. The ratio of V cr n of orientations one and two (denoted by \\ 
and V 2 ) expanded in a power series in 6 is expressed as: 


V 2 

Vi 


1 +2.2 6 + 1.336 2 


(2.4) 


This equation illustrates that small differences in location of the center of mass result in 
'significant variations in impact reorientation velocity. 

To verify the equations and to get a measure of experimental error, tests were per¬ 
formed. Cylinders of different size were placed on a moving conveyor belt. The cylinders 
impacted a steel wall .889 mm (.035 inch) high. Fifteen trials were made for each velocity 
set point and the number of pegs which reoriented was counted. The velocity of the belt 
was monitored by a tachometer which was friction driven directly by the belt. 

Several tests were performed. First, the lower threshold velocity, was determined 

by finding the maximum velocity at which no pegs fell over in fifteen trials. The upper 
threshold velocity, Vf a u, was determined similarly by finding the lowest velocity at which 
all pegs fell over. The median velocity V me d is simply the center of the velocity band 
defined by Vj^jj and V/ a n. Finally, two percentages have been determined from this data. 
The percentage of error is calculated between the theoretical velocity, f / T/ ieory , and the 
median experimental velocity, V me d. The uncertainty band of the experiment is calculated 
from Vj^jj and Vf a u] this band is represented as a percentage of V me d- 

Results confirm that the system model adequately represents the experimental situ¬ 
ation. All tests performed had narrow uncertainty bands. This indicates that there is a 
sharp cutoff between not having enough energy to cause reorientation and having enough 
energy to guarantee reorientation. In addition, the error percentages between Vrheory and 
V mec [ are small. Furthermore, the scaling of the velocity agrees with the scaling predicted 
in equation 2.3. Two pegs which vary greatly in size and length to diameter ratio can have 
the same reorientation velocity. Table 1 gives typical results for two different cylinders. 

The ability to discriminate between orientations, or to reorient one object while not 
reorienting another, depends heavily on the quality of the velocity source. Fortunately, 



excellent velocity sources are inexpensive and easy to build. Figure 2b presents a velocity 
line with the reorientation velocities for the two parts shown in figure 2a. The dark section 
of the line identifies the desired operating velocity range. Within this range part 1 will 
always be reoriented, while part 2 will never be reoriented. By cascading several of these 
impacts at different velocities, a set of like parts, initially oriented randomly can all be 
driven to the same orientation. 


3 Vibratory Reorientation 


3.1 Modeling 

The model for this section is an object in a gravity field. The table on which the 
object rests oscillates vertically (Figure 3). The object is first given some initial angle, 4> o, 
possibly by a sudden motion of the table. For a proper choice of X(f), the table driving 
function, the given object can be maintained in a rocking mode of fixed amplitude. The 
driving function of the table can then be changed, causing an increase in the amplitude at 
which the object will continue rocking. Finally, the amplitude can be increased until the 
object reorients. 

Because of the discontinuity of motion as the part impacts the table, the solution to 
the equations governing the system must be broken into regions. First, the region in which 
there is no impact will be examined. The equations of this motion may be generated using 
Lagrange’s method. The potential energy function can be expressed as 

V = mg(r sin<f) + x) , (3.1) 


where V is the potential energy, m is the part mass, r is the distance from the contact 
point to the center of mass, <j> is the angle between the horizontal and the vector r, and x 
is the vertical position of the table. The kinetic coenergy can be expressed as: 



1 dxd<j> 

H—m — + r— cos <p 

2 dt dt 


1 

H —m 

2 


dcf) 

r—— sin <p 
dt 


(3.2) 


where I c is the moment of inertia of the part about the center of mass. These two terms 
can be used to form the Lagrangian which is then substituted into Lagrange’s equation, 
yielding the two equations of motion for the system, 

, 2 \ d 2 4> , ( , d 2 x 

(I c + mr i ) + mrcos<t> 1 g + 

d 2 x d 2 4> d<f> 2 r 

m —-+ mr—— cos <f> — mr— sin (p + mg — r , 
dt 2 dt 2 dt 

The first of the equations gives the basic motion of the object, while the second equation 
is an expression for F, the contact force. These equations are valid in the conservative 
regions without impact. 
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The next step in formulating the model is to represent the impact of a part with the 
table. At impact, the center of mass of the object has a vertical velocity which is reflected 
with a coefficient of restitution, e. In a perfectly elastic impact, e is equal to one. Physical 
experiments were performed to determine the coefficient of restitution for aluminum parts 
rocking on an aluminum plate. The coefficient of restitution was determined to be between 
.7 and .75 for these tests. For all of the simulations in this paper, an e of .5 was used to 
be conservative. 

As a part hits the table, the impact will also affect the horizontal velocity. If the table 
is frictionless, the part will slide. In this derivation the table was assumed to be rough 
enough so that the contact surface will not slide on the table and after impact the part 
will continue to rock by rotating about the new contact point. 

The next step in creating steady rocking motion was to determine table acceleration 
profiles which would cause a part to oscillate at fixed amplitude. The amplitude of vibration 
is coupled with the frequency of oscillation. As the part gains energy, both the amplitude 
and the period of oscillation become greater. This response makes analysis difficult; few 
tools are available for predicting the behavior of this type of system (Meirovitch, 1975). 

An analytical solution to the system equations for the time required for passage from 
an initial to a final angle can be derived for the case in which the driving function, -ggt 
remains constant. The derivation of this solution is based on conservation of energy in 
the region without impact. For a constant table acceleration, the part can be considered 
to be in a conservative field as shown in equation 3.3a. In this particular case the energy 
balance equation becomes 

1 / o \ d0o / d x . . , 

5( / ‘ + mr )~df +"■(»+ 

= i ± m(g+'^)rsm<t’ , (3.4) 

where 0o is the starting angle of the object, and 0 is some angle of interest. Note that the 
derivative of this expression is the result obtained in equation 3.3a. Equation 3.4 can be 
solved for 


d(f> 

dt 


2m{g + )r 
I c + mr 2 ^ 


(J c +mr 2 )^ 2 
2 m(g+ ^#)r 


+ sin — sin0 


(3.5) 


This equation can in turn be solved by separation of variables and then once integrat¬ 
ing with respect to time. For convenience the constant A is defined as 


A = 


(f,, + 

2 m{g + §f)r 


+ sin 0o 


(3.6) 


and the result is 


f 4 ’" d<t> f T J 2m{g + 

\- A — sin 0 Jo V Ic + mr 2 
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where 4> e is the ending angle of interest. The right side is a simple integral; the left side is an 
elliptic integral. The solution to this equation will be presented as a combination of elliptic 
functions. Though elliptic functions can not be expressed in terms of elementary functions, 
they are well understood and are easily computed. In order to show that equation 3.7 
becomes an elliptic integral, a change of variables must be made. Let 


w = 4> + 


7T 

2 


(3.8) 


so that 


sin 4> = — cos w 


(3.9) 


and 


cos (j)d(f> = sin wdw 


(3.10) 


From these we get 


sin wdw 

d<p = - = dw 

V 1 — cos 2 w 


(3.11) 


Substituting equations 3.11 and 3.8 into the left hand side of equation 3.7 produces the 
following expression 



d<p 

y/A — sin 0 



dw 

\/ A + COS w 


(3.12) 


In the case of impact 


(t>e = i> , 


where ip is defined as the angle between the bottom of the part and the vector r from 
the contact point to the center of mass. If the integral is then rewritten as two integrals 
starting from zero, equation 3.12 can be expressed in terms of the elliptic integral of the 
first kind (Gradshteyn and Ryzhik, 1980): 


l 


4>,+ f 


dw 


y/ A + 


cos w 


Jo 


4 > o + ■?■ 


dw 


VA + cos w 

{ y/2{F{xx,kx) - F(x 2 ,ki)) A< 1 
— {F{x 3 ,k 2 ) - F(x 4 ,k 2 )) A> 1 




. (3.13) 


1 + A 


The elliptic integral of the first kind expressed in normal trigonometric form is 


F(x,k) = [ 
Jo 


arcsm x 


da 


y/l - k 2 sin 2 a 


(3.14) 
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and the limits are 


x\ = 

x 2 = 

ki = 

X3 

X4 

ko = 


1 + sin (j> 0 


A 


1 + sin 4 
1 + A 


1 + A 


4> 0 n 

2 4 

7T 

2 4 


in terms of the system parameters. 

3.2 Choice of Driving Function 


(3.15) 


»2 

As mentioned previously, this analysis is correct only for constant values of The 
first driving function to be considered was a periodic rectangular wave in acceleration (fig¬ 
ure 4). Four parameters must be selected to fully specify the waveform— the acceleration 
upward, AXl; the acceleration downward, AX2; and the switch times, Tl and T2. The 
selection of these parameters is constrained by several relations. First, the areas A1 and 
A2 must be equal in order for the table to return to its initial position and velocity after 
each cycle. Second, the part must start and finish the cycle at the same angle (the sta¬ 
bility requirement). Lastly, the acceleration downward must never exceed gravity. This 
constraint is set so that the part does not lose contact with the table. 

Equations 3.13 and 3.15 were substituted into a nonlinear solution program which 

calculated accelerations and switch times that brought a chosen part back to its initial 

angle in exactly one cycle. Figure 5 shows a time history of the acceleration of the table, 
•2 

, and of the angle <f> as determined by numerical simulation of equation 3.3 including the 
model for impact. For this case, the rocking motion is stable. However, small parameter 
variations cause an instability. Figure 6 shows the rocking motion for an object given a 
slightly larger initial amplitude than that of figure 5. After the first impact, the peak 
amplitude is too low; after the second, it is too high. Eventually, the phasing between 
input and output is lost and the part ceases to rock. 

Because our goal is to bring parts to a stable rocking motion from a range of ini¬ 
tial conditions, the rectangular acceleration waveform will not do the job. A variety of 
waveforms were tested. One which exhibited good performance is shown in figure 7. The 
waveform is simply the rectangular waveform with the edges sloped. Figure 8 shows results 
from a numerical simulation using this forcing function as input. A part is started with 
an initial angle which is twice the steady state angle yet it still reaches a stable equilib¬ 
rium motion. This suggests that by using the proper waveform, a great range of initial 
conditions of the part can be tolerated. 
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If this method is to be used to selectively reorient parts, different parts should exhibit 
significantly different responses. Figure 9 shows variation in response for two similar parts 
given the same initial angle and subjected to the same acceleration profile. For these 
simulations, upward acceleration was set to occur around the time of impact. Several 
simulations were executed with various length “windows” of upward acceleration. The 
ability to discriminate degraded as the window became very large in comparison with the 
difference in impact times. 

An attempt was made to design an experimental apparatus for study of the behavior 
described above. No satisfactory configuration could be created at a reasonable cost. To 
drive reasonably large parts with periods of oscillation greater than .5 seconds would require 
a table with a range of motion substantially greater than that of available electromagnetic 
shakers. Hydraulic shakers are designed to deliver large forces and have servo valves which 
are too small to supply the flow rates necessary to achieve the required velocities. The 
ideal driving system would be a small diameter hydraulic actuator with a relatively large 
servo valve. 

4 C onclusion 

This study of part motion has demonstrated that small variations in inertial properties 
from one part to another cause significant differences in the dynamic behavior of the two 
parts. Two example techniques which capitalize on this property were presented. Practical 
implementation of the first technique, impact reorientation, is feasible given current design 
technology. The second technique, vibratory reorientation, will require construction of 
oscillating tables with performance characteristics available only through custom design. 
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Fig. 1. Impact reorientation model. 
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Fig. 5b. Square wave input. 
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